MgO和Fe2O3的添加对GDCSi体系微观结构及电化学性能的影响

采用溶胶凝胶法制备Gd_0.2Ce_0.8O_3-δ +0.05%（质量分数）SiO₂（GDCSi）电解质。在GDCSi体系中加入Fe₂O₃及MgO可达到降低烧结温度的同时提高晶界电导率，并减小杂质SiO₂对氧离子在晶界处传输的阻碍的目的。将MgO和Fe₂O₃单掺杂或双掺杂在GDCSi体系中并对GDCSi基电解质的微观形貌及电性能进行研究。结果表明，所有样品主要由立方萤石结构相组成；物质的量分数4%MgO单掺杂的GDCSi-M、物质的量分数4%Fe₂O₃单掺杂的GDCSi-F以及物质的量分数2%MgO-物质的量分数2%Fe₂O₃共掺杂的GDCSi-MF均可促进GDCSi体系晶粒增长，降低晶粒间孔隙率，提高电解质的相对密度，降低晶粒电阻R_gi、晶界电阻R_gb及总电阻R_t；GDCSi-MF具有最高晶界电导率和总电导率，在400 ℃时GDCSi-MF的晶界电导率σ_gb和总电导率σ_t分别是GDCSi的10.41和1.82倍。     

关于不定方程x3-1=749y2

利用同余式、平方剩余、Pell方程的解的性质、递归序列证明了:不定方程x³-1=749y²仅有整数解(x,y)=(1,0).     

Moving least squares reconstruction for sharp interface immersed boundary methods

We propose a new approach for reconstructing velocity boundary conditions in sharp-inerface immersed boundary (IB) methods based on the moving least squares (MLS) interpolation method. The MLS is employed to not only reconstruct velocity boundary conditions but also to calculate the pressure and velocity gradients in the vicinity of the immersed body, which are required in fluid structure interaction problems to obtain the force exerted by the fluid on the structure. To extend the method to arbitrarily complex geometries with nonconvex shaped boundaries, the visibility method is combined with the MLS method. The performance of the proposed curvilinear IB MLS (CURVIB-MLS) is demonstrated by systematic grid-refinement studies for two- and three-dimensional tests and compared with the standard CURVIB method employing standard wall-normal interpolation for reconstructing boundary conditions. The test problems are flow in a lid-driven cavity with a sphere, uniform flow over a sphere, flow on a NACA0018 airfoil at incidence, and vortex-induced vibration of an elastically-mounted cylinder. We show that the CURVIB-MLS formulation yields a method that is easier to implement in complex geometries and exhibits higher accuracy and rate of convergence relative to the standard CURVIB method. The MLS approach is also shown to dramatically improve the accuracy of calculating the pressure and viscous forces imparted by the flow on the body and improve the overall accuracy of FSI simulations. Finally, the CURVIB-MLS approach is able to qualitatively capture on relatively coarse grids important features of complex separated flows that the standard CURVIB method is able to capture only on finer grids.     

An improved particle shifting technique for incompressible smoothed particle hydrodynamics methods

The smoothed particle hydrodynamics (SPH) method is one of the powerful Lagrangian tools for modeling free surface flows. However, it suffers from particle disorder, which leads to interpolation and numerical errors. To overcome this problem, several techniques have been introduced until now, among which the particle shifting technique (PST) based on Fick's law is an efficient one. The current form of this method needs tuning parameters to fulfill numerical stability criteria. In this study, to eliminate calibration factors, a new shifting coefficient is derived theoretically based on particle positions before and after shifting, regardless of other parameters such as velocity, pressure, time step intervals, etc. The only required input is particle positions, and the main concern is conserving particle densities in their updated positions. In addition to the proposed PST, a new distribution index (DI) is introduced for measuring the spatial uniformity of particles. Furthering the research, some novel treatments are also studied to improve particle movements near free surface boundary. The proposed idea is only assessed for ISPH method in this study, and its performance in other SPH schemes needs more investigations. Following this innovative method, it is validated by modeling different cases including dam break flow, paddle movement, and elliptical water drop. In all cases, particle arrangements have been improved by means of the modified shifting method. In that sense, good agreements between simulation results with experimental data, analytical solutions, and other numerical methods approve the ability of the developed method in simulating free surface flows.     

A discontinuous Galerkin immersed boundary solver for compressible flows: Adaptive local time stepping for artificial viscosity–based shock-capturing on cut cells

We present a higher-order cut cell immersed boundary method (IBM) for the simulation of high Mach number flows. As a novelty on a cut cell grid, we evaluate an adaptive local time stepping (LTS) scheme in combination with an artificial viscosity–based shock-capturing approach. The cut cell grid is optimized by a nonintrusive cell agglomeration strategy in order to avoid problems with small or ill-shaped cut cells. Our approach is based on a discontinuous Galerkin discretization of the compressible Euler equations, where the immersed boundary is implicitly defined by the zero isocontour of a level set function. In flow configurations with high Mach numbers, a numerical shock-capturing mechanism is crucial in order to prevent unphysical oscillations of the polynomial approximation in the vicinity of shocks. We achieve this by means of a viscous smoothing where the artificial viscosity follows from a modal decay sensor that has been adapted to the IBM. The problem of the severe time step restriction caused by the additional second-order diffusive term and small nonagglomerated cut cells is addressed by using an adaptive LTS algorithm. The robustness, stability, and accuracy of our approach are verified for several common test cases. Moreover, the results show that our approach lowers the computational costs drastically, especially for unsteady IBM problems with complex geometries.     

黏弹性接触界面端附近的奇异应力场

研究蠕变加载条件下线黏弹性材料接触界面端附近的奇异应力场问题.考虑接触界面的摩擦,假设界面端的滑移方向不改变,相对滑移量微小,且其与位移同量级,由此线性化局部边界条件,根据对应原理得到Laplace变换域中的界面端应力场,导出时域中奇异应力场的卷积积分表达式.对卷积积分核函数进行数值反演,考虑接触材料的两类组合,一是持久模量具有量级上的差异,另一是持久模量接近相同.算例结果证实核函数可以用准弹性法求得的解析式较准确地近似.在此基础上,利用积分中值定理,并引入各应力分量的修正系数,得到黏弹性奇异应力场的简化式.结合核函数的数值反演结果分析修正系数表达式的取值范围,得到如下结论,若两相接触材料的持久模量相差很大,可以采用准弹性解的解析式较准确地描述界面端的奇异应力场;一般情况下,应力场不存在统一的奇异值和应力强度系数,当采用类似于准弹性解的表达式近似给出黏弹性应力场时,可以估计此近似描述的误差限.文中最后采用有限元分析黏弹性板端部嵌入部位的应力场,算例包括了黏弹性板与弹性金属支承、黏弹性板与黏弹性垫层所形成的滑移接触界面端,利用黏弹性有限元的数值结果验证理论分析所得结论的有效性.      

正负双向脉冲式爆炸及其诱导的簇发振荡

簇发振荡普遍存在.探索通向簇发振荡的可能路径是簇发研究的热点问题之一."脉冲式爆炸(pulsed-shaped explosion,PSE)"是一种最近被报道的可以诱发簇发振荡的新机制,其特征为平衡点和极限环表现出了与参数变化相关的脉冲式急剧量变.PSE会导致系统轨线急剧跃迁,从而诱发典型的簇发振荡.然而,目前报道的PSE中仅含有"单向的尖峰",未发现"双向的尖峰",且由其诱发的簇发振荡仅含单向的振荡簇.本文以多频激励Rayleigh系统为例,旨在揭示PSE的不同表现形式以及与此相关的簇发动力学.利用频率转换快慢分析法得到了Rayleigh系统的快子系统和慢变量.针对快子系统的分析表明,PSE表现出了较为复杂的动力学特性,其特征是PSE包含了正负双向两个不同的尖峰,此即所谓的正负双向PSE.其急剧量变行为,导致了系统轨线在单个振荡周期内出现正向和负向的多次跃迁,由此得到了由正负双向PSE所诱发的簇发振荡.根据吸引子类型分别揭示了点--点型和环--环型两类簇发振荡模式的产生机制.本文的研究给出了PSE的不同表现形式,丰富了多时间尺度下的簇发振荡的诱发机制.      

基于时域解耦算法的多液舱浮式结构物运动模拟

对于带有多个晃荡液舱的浮式结构物, 浮体的运动、外场水动力以及各舱内的液体晃荡力会实时相互决定, 发生复杂的耦合作用. 为准确模拟多液舱浮式结构物的运动, 本文引入一种有效的时域解耦算法. 该方法以模态分解法为基础, 通过对浮式结构物所受外域水动力和各液舱内非线性晃荡力进行模态分解, 最终形成时域解耦运动方程, 无需迭代求解过程即可显式计算浮式结构物的瞬时加速度. 该方法可避免传统迭代求解方法在迭代次数、截断误差和收敛特性等方面的不足, 减少解耦过程的计算耗时. 本文进一步结合边界元数值方法, 分别对单液舱浮式结构物和多液舱浮式结构物的工况开展数值模拟研究. 通过与单液舱浮式结构物的实验结果对比, 验证了本文时域解耦算法的有效性. 本文详细分析了晃荡力对单液舱浮式结构物运动的影响, 发现存在一个共振影响区间: 当外场波浪频率在该区间之外时, 可以在时域计算结果中观察到稳定的浮体运动; 在比该区间更低频的波况下, 液舱晃荡力与外场波浪力相位相反甚至可以相互抵消, 此时晃荡液舱的存在可以减弱浮体运动; 在比该区间更高频的波况下, 液舱内晃荡力与外场波浪力可以具有相同相位, 此时晃荡液舱的存在会加剧浮体的运动. 本文进一步研究了四液舱浮式结构物在波浪中的纵荡、垂荡和纵摇运动情况, 发现非线性液舱晃荡可对纵荡和纵摇运动产生影响, 但对垂荡运动影响很小.      

人行荷载下梁式结构振动分析的DQ-IQ混合法

为了评估人行荷载作用下梁式结构的振动舒适度，利用微分求积-积分求积，即DQ-IQ混合法求解移动荷载作用下梁的振动响应。人行荷载作用下梁式结构的振动控制方程是含Dirac函数的偏微分方程，首先利用IQ法离散与时间相关的Dirac函数，再利用DQ法把控制方程转化为二阶常系数微分方程，最后利用Newmark算法求解微分方程。以某钢结构连廊为例，利用DQ法计算结构自振频率并与解析解进行对比，结果验证了节点选取和边界条件施加的合理性，再利用DQ-IQ混合法和振型叠加法分别计算了不同行走步频下连廊的响应，计算结果表明，DQ-IQ混合法具有较高的可靠性和精确性。DQ-IQ混合法也可以推广到诸如车辆荷载作用下路面或桥梁的动力响应等其他移动荷载下结构的振动分析。     

不可压粘性流体浸入边界直接力法及软件实现

浸入边界法通过在N-S方程中施加体积力模拟不可滑移固壁边界及动边界，避免生成复杂贴体网格及动网格，极大地节省了网格建模时间及动网格计算消耗。本文提出一种新型附加体积力简化计算方法，将简化附加体积力以源项形式嵌入动量方程迭代中，通过用户自定义函数对CFD软件FLUENT二次开发，实现了浸入边界法和通用流体力学求解器的耦合计算。通过静止圆柱和动圆柱绕流数值模拟进行了验证，并探讨了插值函数对计算精度的影响。研究表明，通过引入浸入边界模型，能够提高计算效率，并实现结构网格背景下复杂边界和动边界的高效建模。